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ABSTRACT 



We confirm discovery by Cao that in the vicinity of fast rotating black holes jets can be launched centrifugally by 
cold, magnetized disks even for nearly vertically shaped magnetic flux surfaces. Outflows produced under such extreme 
conditions are investigated via studying kinematics of test particles in the force-free magnetosphere approximation. 
Implications of a possibility of magneto-centrifugal launching of very well coUimated central outflows around the fast 
rotating black holes are discussed in the general context of the jet formation scenarios in AGNs. 
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1. Introduction 

According to a popular class of astrophysical jet models 
they are powered by rotational energy of accretion disks and 
mass loaded via 'magnetocentrifugal' forces. For Newtonian 
accretion disks such models require the poloidal magnetic 
field lines to be inclined to the geometrically thin accre- 
tion disk by less than 60 degrees, independently on a dis- 
tance from the central object (Blandford & Payne 1982). 
For disks around black holes (BHs) this angle depends on 
the radius and the BH spin (Cao 1997; Lyutikov 2009). It 
is the larger the smaller is the distance and the larger is the 
BH spin. At the inner edge of a Keplerian disk (located at 
the marginally stable orbit - r,„s) and for the maximal BH 
spin it approaches 90°. We confirm these results and inves- 
tigate the kinematics of test particles in the fixed force-free 
magnetosphere. Our paper is organized as follows. General 
equations describing the particle kinematics are derived in 
§2. Geometry of the effective equipotential surfaces and of 
"light cylinders" are presented in §3. Kinematics of test par- 
ticles and its dependence on the BH spin and the launch- 
ing distance is illustrated in §4. Results are discussed in a 
general context of production of relativistic jets in active 
galactic nuclei in §5 and summarised in §6. 



2. A rigid rotation of particle trajectories in the 
Kerr metric 

Kinematics of particles forced by magnetic fields to move 
on rigidly rotating trajectories is convenient to investigate 
in a frame co-rotating with magnetic field lines. In such a 
frame the norm of a 4- velocity, m^u' = —1, is in the Kerr 
metric given by equation: 

{gu+^gt^^Q+g^^nDiu')"^ +2{gt4>+g^^no)u*u'''' + (1) 
g^^{u'^')^+grr{un^+gee{uy = -1, 

where (f)' = 4> — riot, u'^ = u'^ — fio^*, = const, the t, 
0, r, and 9 are the Boyer-Lindquist coordinates, and gik = 



gik{fid]a) are the Kerr metric components. Through the 
paper the following quantities are expressed in dimension- 
less units: radius r : r/{GM/c^), BH spin a* = J/{GM'^/c), 
and angular velocity fip : fio/(c'^/GM), where M and J are 
the BH mass and angular momentum, respectively. 

For a particle trajectory which in the rotating frame is 
described by functions r = r{9) and 0' = (f>'{0), the r and 
(/>' components of the particle 4-velocity are = r^gu^ and 

u'^ = 4>' gu^ , respectively. Then the norm of the 4-velocity 
reads 



gtt{u) +2gteuu + gee[u ) =-1, 
where 

gtt — gtt + 2(7(0^^0 + g<}>4>^o , 

gte = [gt^ + g<t><t>^Q)(t>\e ^ 

gee = gee + grr{r,eY + 9<l><l>i(t>' ,ey ■ 



Noting that 

Uf = gttu*- + gteu'^ 



(2) 

(3) 
(4) 

(5) 
(6) 



where e is the constant of motion of particle moving along 
a rigidly rotating trajectory, one can find inserting ut from 
Eq. (6) into Eq. (2) that 



{-gttjgee + {gtef 



and 




(7) 



(8) 



{^Y " (~g») 

{-gtt)gee + {gteY 



Hence, for a fixed particle trajectory and a given con- 
stant of motion e kinematics of a test particle is fully de- 
termined. This kinematics can be illustrated in the locally 
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non-rotating frame (Bardeen, Press, & Teukolsky 1972), i.e. 
the frame of zero-angular- momentum-observers (ZAMO). 
In such a frame the hue element is 

ds'^ = {gtt+i^gt4,)dt'^+g^^{d4}-ujdt)'^+grrd7^^+geed9'^ ,(9) 

where uj = —gt<p/g4>4> is the angular velocity of 'the dragged 
inertial frames'. Then the projection of the 4- velocity onto 
the orthonormal tetrad of the local Minkowski space can 
be used to calculate the Lorentz factor and velocity com- 
ponents of a particle in the ZAMO frame: 



7 = u^*^ = \/-9tt + 9t4,/9<t'<t> ' 



(10) 



{^lo - uj) 



9<P<P 



-9tt + 9t4,/9H 



.(11) 



2 U" / grrTfi + 906 



" V -9tt + 



(12) 



where the physical velocities w*^*-' are expressed in the speed 
of light units. 

3. Effective potential 

An effective potential defined as the minimum energy of test 
particles forced to rotate with a given angular velocity f^o 
is Veff = \f—gtt (obtained from Eq. (7) setting = 0). Its 
equipotential surfaces, T4//(r, 9) = const, are illustrated in 
Fig.[TJ They are enclosed between the inner and outer 'light 
cylinders' given by gu — (Lyutikov 2009). Locations of 
the inner and outer light cylinders in the equatorial plane 
are presented in Fig. [5J For a given spin the cylinders co- 
incide at the photon orbit. For a* > 0.91 and the angular 
velocity fig corresponding to the marginally stable orbit 
the outer light cylinder is enclosed by the BH ergo- 
sphere. In Fig. [3] we show the dependence of the equatorial 
and asymptotic locations of the outer light cylinder on the 
BH spin for fig calculated at the marginally stable orbit. 
The asymptotic radius is always close to its equatorial plane 
value (e.g. 1A.7M vs 13.6Af for a non-rotating BH). 

Equipotential surfaces, dVeff — 0, intersect the equato- 
rial plane at angles 



tan^ 



d0\ 



eff.r 



=7r/2 \^e//,8/ e=7r/2 

Since V;//,e|e=7r/2 = 0, while Veff,r\e=-n/2 = only at r = 
■To at which the Keplerian law is satisfied, i.e. 

^0 - • (14) 

this angle is ^ = 7r/2 for all radii but r — tq. Kt r ^ tq this 
angle can be found by applying in Eq. (13) the L'Hospital's 
rule and by noting that Ve//,re|e=7T/2 = 0- This gives 



(13) 



tan^o 




(15) 



7r/2,r— ro 



1 - 4ar, 



-3/2 







This formula fully agrees with the formula obtained by Cao 
(1997) and Lyutikov (2009). Dependence of on the radius 
ro and on the BH spin is shown in Fig. 01 




4 6 

r sin e [M] 

(a) a* = 
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r sin [M] 
(b) a* = 0.99 

Fig. 1. Equipotential surfaces of the effective potential 
(^e// = V 9tt) for fio ^ floir ^ rms,a*) for a non-rotating 
(top) and spinning (bottom panel) BHs. The thick solid 
lines represent Veff(r,9) = Vq = Veff{r = r,ns,9 = tt/2) 
being the effective potential crossing the equator at r = 
rms- The dashed lines present locations of the inner and 
outer light cylinders. The thin solid lines are drawn for the 
following values of the effective potential: |, |, |, ^Vq. 



4. Kinematics of test particles in a force-free 
magnetosphere 

4.1. Purely poloidal magnetic fields 

For a strong dynamical domination of large scale magnetic 
fields driven by an accretion disk, test particles are re- 
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Fig. 2. Locations of the inner and outer light cyhnders for 
^ = 7r/2 as a function of radius tq defining the angular 
velocity flo- Inner and outer cylinders are represented by 
curves below and above the dashed line r = rg , respectively. 
Dotted curves denote locations of light cylinders for ro < 
Tms- Profiles for three values of BH spin: -0.999, 0, 0.999 
are presented. 
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Fig. 3. Locations of the light cylinders for Oq = ^air — 
fms) versus BH spin. Location of the outer light cylinder is 
shown for the equatorial plane (thick, dashed line) as well 
as its asymptotic {0 0) value (thin, dashed line). 
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Fig. 4. The critical angle dependence on the radius rg. 
Profiles for different BH spins, starting at the marginally 
stable orbits, are shown. 
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Fig. 5. Dynamics of a particle moving along a purely 
poloidal magnetic field lines achored at r^s and inclined 
at different angles for non-spinning (solid) and spinning 
(dashed lines) BHs. The horizontal axis represents the ver- 
tical coordinate z. 



stricted to move along magnetic field lines. Being anchored 
to a differentially rotating accretion disk, the poloidal mag- 
netic field generates differentially rotating magnetosphere. 
Trajectories of test particles can then be identified with 
magnetic field lines rotating with an angular velocity of 
their foot-points (Ferraro 1937). For purely poloidal mag- 
netic fields such trajectories are in the co-rotating frame 
planar (0' = const), the metric component gte = and 
equations (7) and (8) reduce to 



u = 



' (go)^-(-gtt) 

{-gtt)999 



Co 



i-9tt) 



(16) 



(17) 



where eg = -gttiO = 7r/2,r = ro), ff^fc = g^k{O,r{0;ro)), 
and r{9;ro) is the shape of the magnetic field surface which 
is determined by the shape of poloidal magnetic field lines. 

Test particles can be pulled from the cold disk by cen- 
trifugal forces provided the inclination of magnetic field sur- 
faces at the foot-point is smaller than and they can pro- 
ceed further provided the shape of poloidal magnetic field 
lines is such that along them the condition —gu < ^ is 
satisfied. 

Examples of kinematics of the test particles are plotted 
in Fig. 4. The illustrated cases are calculated for magnetic 
field lines anchored at the marginally stable orbit and hav- 
ing shape 



r = ro 



tan^o 



(sin 6* tan ^ — cos 6*) ' 



(18) 



which in the coordinate plane (r cos 9, r sin 9) is the straight 
line anchored at r = ro and inclined at the angle ^. The 
presented velocities and Lorentz factors were calculated in 
the locally non-rotating frame (see Eqs. (10)-(12)). As we 
can see from Eq. (8), the particles which are forced to move 
on rigidly rotating planar trajectories approach the speed 
of light at the light cylinder. This simply demonstrates that 
no physical solutions exists on and beyond light cylinders 
for purely poloidal magnetic field structures. 
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4.2. Force-free magnetosphere with non-zero azimuthal 
magnetic field component 

Not necessarily such constraints apply for magnetosphere 
if taking into account presence of toroidal magnetic field 
component which is induced by poloidal electrical currents. 
In this case the Ferraro's law about iso-rotation still ap- 
plies, i.e. magnetic field lines and, therefore, particle tra- 
jectories, still rotate rigidly, however, due to the presence 
of the toroidal magnetic component the rotating trajecto- 
ries are not planar and the test particles may pass the light 
cylinder sliding on magnetic field lines in the opposite di- 
rection to the rotation of magnetic field surfaces, i.e. with 
d(f>' /dt < 0, which takes place if 0' g < 0. This can be con- 
cluded from Eq. (8) after rewriting it in the form 



(-5tt) 



1 - 



' i-(-g»)A1 
1 + {~-gtt)g6s/9t0 



(19) 



where the metric functions gik are calculated along trajec- 
tories r = r{d;ro) and (j>' — (p'{d;rQ). At the light cylinder 
both denominator and numerator {—gtt and the expresion 
in brackets) become zero and applying the L'Hospital's rule 
one can check that u* and, therefore, 7 given by Eq. (10) 
is finite. Hence, the test particle can pass the light cylinder 
and proceed further provided 



(ffte)^ > gttgee , 

i.e. for respectively large negative values of </)' , 



gttjgee + grrjrfiY) 

{gt<p + g<i,^^oY - g<i,^gtt 



(20) 



(21) 



Examples of kinematics of test particles for helical mag- 
netic field structure are illustrated in Fig. [B] The veloc- 
ity components and the Lorentz factor are computed for 
r{9;ro) given by Eq. (18) and for 



cos^ 



(r sin 6 — tq) 



(22) 



The above choice of (/)'(6';ro) corresponds roughly with 
a poloidal magnetic field component scaling as Bp ^ 
l/(rsin6')^ and toroidal component scaling as ^ 
l/{rsmO), while 77 w -Bp/ B^\r sin e=i /no (in Fig. [7] we plot 
a visualisation of such a shape of the magnetic field lines) . 

The upper panel of Fig. |6] presents the velocity and 
Lorentz factor profiles for a non-spinning BH. The par- 
ticle is initially at rest with respect to the disk — its 
poloidal velocity component in the ZAMO frame is zero 
while its azimuthal velocity corresponds to the Keplerian 
angular velocity at tq. Once the particle leaves the equilib- 
rium it starts to be centrifugally accelerated along the mag- 
netic field line. The poloidal velocity component rapidly 
increases. Due to the fact that the magnetic field lines are 
inclined also in the azimuthal direction the particle starts 
to slide along them in the direction opposite to the disk ro- 
tation. Therefore, its angular velocity decreases below the 
initial value fig. The corresponding profile of the physical 
velocity in the azimuthal direction shown. The light 

cylinder is crossed at z « 14M with both velocity compo- 
nents close to 0.5c. The profile of the Lorentz 7 factor is 
shown with the solid line (please note the vertical scale for 
7 is marked on the right axis). Initially, it corresponds to 
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Fig. 6. The velocity components and the Lorentz factor 
7 of particles moving along magnetic field lines with non- 
zero azimuthal component anchored at Vms and inclined at 
the critical angle ^ = ^0 for non-spinning (top) and spin- 
ning (bottom panel) BHs. Poloidal and azimuthal velocity 
components are marked by dotted and dashed lines, respec- 
tively. Solid curves present the Lorentz factor profiles. 



the Keplerian angular velocity ar r^s (~ 1.2 for a* = 
case). At the light cylinder 7 w 1.7. The bottom panel of 
Fig. [5] presents a similar study of particle dynamics for a 
rotating BH case. Its behaviour is similar with exception to 
the fact that the 7 factor is initially strongly dominated by 
the azimuthal velocity component. 

We demonstrate in Fig.[5]the dependence of the test par- 
ticle kinematics on ^, tq and rj. The upper panel shows that 
the lower the inclination angle ^ the more rapid is the parti- 
cle acceleration. The middle panel presents the dependence 
on the r] parameter which determines the magnetic field 
line torsion. The more 'twisted' is the line (corresponding 
to increasing rj) the easier particle slides in the azimuthal 
direction and the resulting 7 factor increases slower. Finally, 
the lower panel presents the dependence on the radius 
defining the location of the magnetic field line foot-point. 

Obviously, illustrated in this section kinematics has only 
exemplary character and applies to the region in which 
the structure of the outfiow is strongly dominated by the 
magneto-dynamics. For low mass loading rates, such a 
force-free approximation can be used to study the matter 
kinematics even far beyond the light cylinder. However, for 
high mass loading rates the force-free approximation can 
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Fig. 7. Shape of the magnetic field hues with rj — 1.0 for 
^ = 30° (the thickest sohd hue), ^ = 45° and ^ = 60° 
(the thinnest sohd hne). For the latter (^ = 60°) two other 
curves are presented for r] = 0.95 (dotted line) and rj = 1.05 
(dashed line). All lines are anchored at r^s and the BH 
has zero spin. The black spot denotes location of the BH 
horizon. 

break down even in the sub-Alfvenic region, deep within 
the light cylinder. 

5. Discussion 

Magnetocentrifugal scenario is applicable to the production 
of jets in 'Newtonian' objects like YSOs (Young Stellar 
Objects). However, it cannot explain the whole jet pro- 
duction process in BH accretion systems. This is because 
at least in idealized, steady, axisymmetric models with an 
accretion disk treated only as the boundary surface, effi- 
cient centrifugal mass loading makes the outflows too heavy 
to be accelerated up to the observed relativistic velocities. 
Favored mass loading scenarios for such objects involve pair 
production by photon-photon interactions. For typical BH 
accretion systems the coronal activities are enough power- 
ful to support in the vicinity of the BH density of pairs by 
many orders of magnitude larger than the Goldreich- Julian 
density (Phinney 1983; McKinney 2005), which allows to 
treat the outflows in the ideal MHD approximation. By the 
same time, the rest mass-energy density of the pair plasma 
is by several orders lower than energy density of magnetic 
flelds required to power luminous large scale jets in radio- 
loud objects, and this implies formation of relativistic, at 
least initially strongly Foynting flux dominated, outflows. 

The e"'"e~ jets can be powered by both rotating BHs 
(Blandford & Znajek 1977; Phinney 1983; Beskin 1997) and 
accretion disks (Blandford 1976; Lovelace 1976; Lovelace, 
Wang & Sulkanen 1987; Ustyugova et al. 2000). Such jets 
need external collimation to be effectively accelerated and 
converted to the matter dominated ones (Begelman & Li 
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Fig. 8. The Lorentz 7 factor dependence on the magnetic 
fleld line inclination angle ^ (top panel), rj parameter (mid- 
dle panel) and the vq radius (bottom panel) for a non- 
spinning BH. The vertical dashed lines denote crossing 
through the outer light cylinder. If not stated otherwise 
the parameters are: ^ = 60°, 77 = 1.0 and r = Vms- 

1994; Narayan , McKinney, & Farmer 2007; Komissarov, 
et al. 2007; Lyubarsky 2009; Forth & Fendt 2010). The 
collimation is often considered to be provided by slower, 
barionic MHD outflows launched in the accretion disk 
(Sol, Pelletier, & Asseo 1989; Bogovalov & Tsinganos 2005; 
Gracia et al. 2005; Beskin & Nokhrina 2006). However, such 
a jet structure may need modification for jets launched 
around fast rotating BHs if the possibility of centrifugal 
launching of proton-electron outflows from innermost por- 
tions of accretion disk was taken into account. Since the 
inclination angle of the effective potential surfaces is max- 
imal at the inner edge of a disk and rapidly drops with a 
radius (see Eq. 15 and Fig. |4l the p-e dominated outflow 
from the innermost portions of accretion disk is embraced 
by the Foynting flux dominated electron-positron outflow. 
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This implies less relativistic spines and more relativistic 
sheaths of jets produced by the accretion disk around fast 
rotating BHs. or a respective triple jet structure if taking 
into account the central contribution to the outflow from 
the BH magnetosphere. 

The production of powerful jets by accretion disks, both 
centrifugally loaded and by pairs, requires strong large scale 
magnetic fields. Possibility of dragging of such fields by ac- 
creting plasma, originally suggested by Bisnovatyi-Kogan & 
Ruzmaikin (1974), is still debated (Lubow, Papaloizou, & 
Pringle 1994; Spruit & Uzdensky 2005; Bisnovatyi-Kogan 
& Lovelace 2007; Rothstein & Lovelace 2008; Beckwith, 
Hawley, & Krolik 2009). Related to the lack of the pro- 
found model of the evolution of magnetic fields in accretion 
disks, unknown remains radial distribution of inclination of 
magnetic field lines which is critical to establish a launch- 
ing distance domain of centrifugal outflows and their initial 
colhmation. And finally, even for a fixed large scale mag- 
netic field intensity and geometry, mass loading rate and 
therefore the terminal speed of the centrifugal outflow de- 
pends very much on details of the vertical disk structure 
(Ogilvie & Livio 2001), which due to severe uncertainties is 
usually ignored. In particular, due to partial losses of angu- 
lar momentum (taken away by the outflow) the boundary 
layers of the accretion disk likely become sub-Keplerian. 
As a result, proton-electron outflow may need some initial 
boost to pass the created potential barrier. It can be pro- 
vided by heating or mechanically by some magnetic flar- 
ing activities and/or by radiation pressure of effectively 
super-Eddington flux. Existence of the additional potential 
barrier may significantly limit mass loading rate allowing 
proton-electron loaded outfiows to reach at least mildly rel- 
ativistic speeds. The resulting jet structure — mildly rela- 
tivistic proton-electron component sandwiched between the 
pair dominated relativistic spine and sheath — albeit very 
speculative at the moment, is very promising from the ob- 
servational point of view because may explain a significant 
proton content of AGN jets deduced from analyses of a mat- 
ter content in blazars (Sikora & Madejski 2000; Ghisellini 
et al. 2009) and in radio lobes of powerful radio galaxies 
(Stawarz et al. 2007; Pcrlman et al. 2009). 



6. Conclusions 

Main results of this paper can be summarized as follows: 

- An effective potential and light cylinders are investigated 
in the rigidly rotating frame in the Kerr metric. The in- 
tersection of equipotential surfaces with a geometrically 
thin disk at the annulus where a given angular velocity is 
equal to the Keplerian velocity gives the critical angle below 
which a cold outflow can be launched by centrifugal forces. 
The location of outer light cylinders is shown to depend 
strongly on radius and BH spin. For BH spin a* > 0.91 
and the Keplerian angular velocity at the marginally sta- 
ble orbit the light-cylinder radius in the equatorial plane is 
enclosed within the erghosphere. 

- The condition for magnetocentrifugal launching of jets 
obtained for Keplerian disks rotating around Kerr BHs by 
Cao (1997) and Lyutikov (2009) is conflrmed. It tells us 
that the maximum inclination angle of magnetic flux sur- 
faces at which cold matter can be extracted from the disk 
depends on a distance from the BH and on the BH spin. 
The condition shows that in case of very fast rotating BHs 



the central outflows can be launched even along of almost 
vertically shaped magnetic surfaces. 

- We show how kinematics of test particles pulled by cen- 
trifugal forces from a Keplerian disk can be algebraically 
determined for a given magnetic field structure in the force- 
free outflow approximation. Examples of test particle kine- 
matics are illustrated and the condition for the toroidal 
magnetic field component is derived to allow the particle 
to cross the light cylinder. 

- Possible implications for a jet structure are discussed as 
imposed by the condition for magnetocentrifugal launching 
of jets by inner portions of magnetized disks around fast 
rotating BHs. In this case a triple-component structure of 
a jet can be envisaged, with a proton-electron component 
of a jet being sandwiched between the relativistic pair dom- 
inated spine and sheath. 

Acknowledgements. AS and MS acknowledge direct support by Polish 
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